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Analysis  of  propagation  loss  is  generally  based  on  the  absorption  formula 
developed  by  Thorp  in  1965.  The  Global  Model  is  an  improved  formula  based 
on  more  recent  field  and  laboratory  experiments.  The  new  formula  involves 
the  pH  factor  which  varies  with  both  region  and  depth.  In  order  to 

account  for  depth  variation,  the  effective  absorption  must  be  calculated  by 
integrating  loss  over  ray  paths.  Contour  charts  are  provided  for  estimating 
the  required  K  profiles  throughout  the  Worjd  Ocean. 

^  fhe'purpose  df^his  report  ts-td  examin&€he  impact  of  the  model  on  sonar 
performance  and  t'he  surface  and  bottom-loss  models  upon  which  predictions 
are  currently  based.  The  present  analysis  is  limited  mainly  to  active  systems 
operating  in  the  frequency  range  3-4  kHz.  The  four  regions  examined  are: 
North  Atlantic,  North  Pacific,  Eastern  Mediterranean  and  Norwegian  Sea. 

Comparison  of  K  model  predictions  with  those  based  on  the  Thorp  fotmiula 
indicates  that  two-way  path  absorption  for  the  first  CZ  zone  ranges  from 
approximately  4  dB  less  in  the  North  Pacific  to  1 2  dB  greater  in  the  Eastern 
Mediterranean.  In  the  bottom-bounce  mode,  these  values  are  reduced  roughly 
in  proportion  to  range. 

The  MGS  bottom-loss  model  is  based  on  analysis  of  one-way  paths  using 
the  Thorp  formula.  The  difference  between  the  K  model  and  Thorp  formula  for 
I  O'*  grazing  angle  in  the  Eastern  Mediterranean  is  roughly  4  dB,  making  the 
actual  bottom-loss  smaller  by  about  one  province  number.  The  differences 
decrease  with  increasing  grazing  angle  and  are  smaller  in  the  other  areas.. 
Corrections  to  the  MGS  model  should  therefore  be  relatively  minor.  ^ 

K  model  predictions  of  two-way  loss  for  the  surface-duct  mode  at  25  km 
range  vary  from  3  to  8  dB  more  than  Thorp  values,  since  pH  tends  to  be  high 
near  the  surface.  The  higher  absorption  coefficients  evidently  account  for  a 
large  fraction  of  the  excess  attenuation  previously  ascribed  to  surface  loss. 
Since  the  surface-loss  models  were  based  on  pre-Thorp  absorption  models,  a 
major  reassessment  is  indicated 
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Sound  absorption  in  sea  water  is  more  than  an  order  of  magnitude  greater 
than  in  fresh  water.  The  excess  is  the  result  of  relaxations  involving  certain 
chemical  constituents.  Laboratory  experiments  in  the  1950’s  identified  the 
principle  mechanism  as  a  magnesium  sulfate  with  relaxation  frequency  near 
100  kHz.  Combining  measured  relaxation  parameters  with  data  from  field 
experiments  in  the  frequency  range  2-30  kHz,  Schulkin  and  Marsh  (1)  devised 
the  first  practical  formula  for  sound  absorption  in  sea  water  in  1962. 

At  about  that  time,  interest  shifted  to  lower  frequencies  in  an  effort  to 
increase  range  and  questions  arose  about  extrapolation  of  the  S&M  formula. 

In  1965,  Thorp  [2]  reported  results  of  long-range  sound-channel  experiments 
covering  lower  frequencies  and  a  second  anomaly  was  observed.  The  data 
were  fitted  by  adding  a  1  kHz  relaxation  component  to  the  S&M  formula  and 
this  became  known  as  the  Thorp  formula"  (3). 

Mediterranean  experiments  by  Leroy  (4),  showed  a  similar  anomaly  with 
greater  magnitude  and  relaxation  frequency.  Experiments  carried  over  the 
next  two  decades  in  other  areas  confirmed  the  variability  of  the  anomaly  and 
regional  dependence  became  a  critical  factor 

The  anomaly  was  identified  as  a  second  relaxation  involving  boric  acid,  a 
minor  sea-water  constituent  (51.  In  this  case,  absorption  depends  on  pH  as 
well  as  temperature  (61.  The  relaxation  mechanism  was  identified  as  the 
boric  acid/carbonate  equilibrium  and  parameters  were  measured  in  the 
laboratory  171.  The  experiments  also  showed  that  a  pH-dependent  relaxation 
of  magnesium  carbonate  has  a  minor  but  significant  role  in  sea  water 
absorption  (01. 

An  absorption  formula,  based  solely  on  known  chemical  processes,  would 
be  too  complex  and  the  accuracy  would  be  inadequate  as  well.  However,  the 
range  of  the  ocean  parameters  is  small  and  simplifying  approximations  can 
be  made,  i  e  that  magnitudes  of  pH-dependmt  terms  increase  exponentially 
with  pH  and  that  relaxation  frequency  increases  exponentially  with  the 
temperature  The  Thorp  formula  con  then  be  modified  by  simply  adding  the 
third  relaxation  and  including  pH  and  temperature  factors  for  all  terms. 

A  3-relax8tion  formula  has  been  developed  and  absorption  coefficients 
based  on  archival  pH  values  have  been  found  to  be  in  good  agreement  with 
experimental  data  For  prediction  purposes,  the  Global  Model  (91  includes 
contour  charts  for  estimating  pH  correction  factors  vs  depth  as  well  as 
region  The  purpose  of  this  report  is  to  examine  the  potential  impact  of  the 
model  on  systems  performance  and  the  surface  and  bottom-loss  models  upon 
which  predictions  are  currently  based 
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Absorption  Model 


A=A,(MgS04)+A2(B(OH)j)+A3(MgC03) 

A„={S/55)a„F^F„/(F|f+F^) 


flj  -  0.5  XIO 


D(km)/20 


02"  O.lxJO 


(pH- 8) 


83=  0.03xt0 


(pH-8) 


Fj  =  50x10 
F2=  0  9x10^ 
F3  =  4.5x10^ 


T/60 


T/30 


Atlantic  4*c  pH  8.0 

A=0.007f 2  +0  1  f 2/(1  ♦  f 2)+o.i8  *) 


N. Pacific  4*C  pH  7.7 

A  =  0  007f2  ♦0.05f2/(  !♦  f 2)*o  09  f2/(6^*f^) 


Mediterranean  Sea  M*C  pH8.3 
A=0.006f2^0  26f2/(1.4^»f^)^0.7ef^/(12^^f^) 


Norwegian  Sea  0*C  pH  8.2 

A=0.01f*>0.16  f^/(0.9^f^)^0.2  f^/(4.5^*f^  ) 


Table  1;  Simplified  absorption  formulae. 


Relaxational  absorption  formula  nave  the  simple  form  shown  in  the  top  box 
of  Table  I  Values  are  for  A  in  dB/km,  F  and  relaxation  frequency  F^  in  kHz 

and  T  in  *C  The  3  components  of  the  model  ere  simply  additive.  Thorp's  value 
pH=0  0  IS  used  as  reference.  The  pure  water  term  is  omitted,  making  the 
formula  valid  for  frequencies  less  than  roughly  100  kHz.  Values  of  pH  in  the 
World  Ocean  vary  roughly  from  7.7  to  6.3,  coiresponding  'n  absorption 
ratio  of  almost  4/1  at  the  lower  frequencies 
The  magnesium  sulfate  term  (Al)  also  includes  the  depth  factor  D(km), 
which  is  adapted  from  the  pressure  correction  of  Fisher  and  Simmons  (10), 
Depth  dependencies  of  the  other  relaxations  are  not  yet  known,  however,  the 
field  data  indicate  that  effects  on  the  boric  acid  term  (A2)  are  negligible 
Effects  on  the  magnesium  carbonate  term  (A3)  will  be  much  less  important 
because  its  contribution  is  so  small 

Specific  foiTnulae  are  shown  in  the  bottom  box  Note  that  the  Al  terms  are 
approximations  valid  only  for  F«  F,. 
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Figure  1:  Thorp  s  data  and  3-relaxation  model. 

Thorps  formula  for  the  N  Atlantic  sound  channel  was  given  as  (3j: 

Ck:=40  FV(4100>F2)*0.j  f2/(Uf2)  dB/Xyd. 

Note  that  the  temperature  and  depth  terms  in  the  S&M  formula  are  omitted 
and  therefore  it  applies  strictly  to  roughly  4‘’C  at  I  km  depth 
Figure  1  shows  Thorp's  original  data  compared  to  the  3-relax8tion  model 
The  individual  components  are  identified  end  the  top  curve  is  their  sum  The 
overall  fit  to  the  data  is  as  good  or  better  than  the  2-component  formula. 
Note  that  the  boric  acid  (A2)  coefficient  is  lower  than  that  in  the  Thorp 
formula  by  some  lOS,  making  the  coefficient  in  dB/km  becomes  equal  to  the 
coefficient  of  the  second  term  in  Thorp's  equation.  The  magnesium  carbonate 
component  makes  up  the  difference  at  the  lower  frequencies 
The  parameter  adjustment  in  the  present  model  is  mainly  justified  on  the 
basis  of  data-fit,  the  third  component  having  been  found  essential.  When  the 
laboratory  date  were  fitted  by  a  2-rel8xation  model,  serious  discrepencies 
appeared  at  higher  pH  values  For  example,  correcting  Thorp's  second  term  to 
pH0  5  gave  values  that  much  too  low  in  the  10  kHz  range,  clear  evidence  of 
another  component  By  sea-water  synthesis,  the  mechanism  was  identified 
as  magnesium  carbonate,  temperature  and  pH  dependence  of  the  parameters 
were  measured  and  the  relaxation  incorporated  as  the  third  component  of  the 
absorption  model 
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Figure  2;  Model  and  data  comparison 

Figure  2  compares  the  3-component  relaxation  model  predictions  and  data 
from  sound-channel  experiments  in  the  Mediterranean.  North  Pacific  and 
North  Atlantic  (ThorpV 

-he  N:.rth  Pacific  case,  the  lower  value  pH^?,?  reduces  both  the  bone 
acid  (A2)  8i:1  the  magnesium  carbonate  (A3)  coefficients  by  a  factor  of  two 
compared  to  the  N  Atlantic  values  Relaxation  frequency  depends  only  on  the 
temperature  and  remains  the  same 

In  the  Mediterranean  case,  the  higher  value  3  increases  both  the  boric 
acid  (A2)  and  the  magnesium  carbonate  (A3)  coefficients  by  a  factor  of  two 
compared  to  the  N  Atlantic  However,  the  relaxation  frequency  is  higher  and 
the  curves  do  not  differ  by  as  large  a  large  factor  at  the  lower  freq  uencies 

The  value  pH  6  0  has  been  assumed  for  Thorps  sound-channel  experiment 
and  is  used  as  reference  value  for  the  3-relaxation  model  Predictions  based 
on  archival  pH  values  then  show  good  agreement  within  experlmentaf  limits 
for  all  regions  examined 
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Frequency  (kHz) 

Figure  3:  Absorption  spectra  A/F  vs  pH  forT=4®C. 

The  effects  of  the  3-rel8X8tion  model  ere  illustrated  in  Figures  3  and  4. 
Note  that  the  absorption  coefficient  A  has  been  divided  by  the  frequency 
F(kHz)  in  order  to  compress  the  vertical  scale. 

Predicted  spectra  vs  pH  for  4®C  and  normal  salinity  (5-35)  are  shown  In 
Figure  3  and  cover  the  approximate  sea-water  range  of  pH  7.7-8  3  In  0. 1 
steps  The  effect  of  pH  on  the  absorption  spectrum  below  10  kHz  is  mainly 
due  to  the  coefficient  of  the  boiic  acid  component  A2.  Below  1  kHz.  the 
max/min  ratio  approaches  the  factor  four. 

From  analysis  of  experimental  data,  the  overall  error  of  the  model  is 
estimated  to  be  roughly  of  the  order  of  ♦  15S.  In  other  words,  the  RMS  error 
of  an  estimated  coefficient  In  dB/km  is  not  expected  to  be  any  g^^eater  if  not 
limited  by  the  accuracy  of  the  environmental  factors.  For  example,  a  change 
of  0  05  pH  units  corresponds  to  a  ^2%  change  in  the  boric  acid  coefficient 
and  the  pH  error  must  be  within  these  limits  in  order  to  realize  the  model 
accuracy  A  major  limiting  factor  is  therefore  the  accuracy  of  of  the  pH 
data-base. 
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Frequency (kHz) 

Figure  4:  Absorption  spectra  A/F  vs  temperature  for  pHB.O, 

Figure  4  shows  the  predicted  cpectrum  dependence  on  temperature  for 
pH=6  0  and  normal  salinity  (S=35).  covering  the  nominal  sea-water  range 
C-30“C  in  5X  steps. 

Temperature  effects  in  the  model  come  about  only  through  the  relaxation 
frequency  coefficients  of  the  three  terms  Increasing  temperature  shifts  the 
curves  upward  in  frequency 

Model  errors  involving  temperature  coefficient  are  estimated  to  be  within 
the  expected  il5S  limits  Since  the  temperature  sensitivity  is  rather  low, 
the  value  need  only  be  known  only  to  within  roughly  ♦  I®C  to  be  well  within 
these  limits 

Errors  involved  in  the  salinity  factor  ere  probably  also  within  the  accuracy 
limits  for  the  normal  sea-water  range  Salinity  variations  are  evidently  the 
result  of  evaporation  or  dilution  by  fresh  water  Concentrations  of  all  the 
constituents,  except  carbon  dioxide,  tend  to  remain  in  constant  ratio  Carbon 
dioxide  conc6ntrat.cn  evidently  tends  to  remain  constant. 

The  multiplying  factor  S/35  for  all  three  components  may  be  justified  over 
a  limited  range,  say  30-40  ppt  Beyond  this  range,  corrections  must  be  made 
for  changes  in  the  boric  and  retixation  frequency  This  can  require  a  more 
accurate  knowledge  of  the  concentrations,  which  is  beyond  the  scope  of  the 
present  model 


PH  Profiles 


Figure  5  pH  profiles. 

In  addition  to  regional  vanapility,  pH  also  can  depend  strongly  on  depth 
Figure  5  compares  pH  profiles  for  various  regions  of  the  World  Ocean  It  is 
clear  that  the  absorption  will  vary  significantly  for  different  ray  paths, 
particularly  in  the  North  Pacific  where  sound-channel  losses  will  be  much 
less  than  for  other  paths  in  order  to  account  for  the  depth  variability,  the 
most  general  method  of  estimating  effective  adsorption  is  to  integrate  loss 
over  all  the  ray  paths  involved  This  requires  an  appropriate  profile  for  the 
pH  factor  The  Global  Model  uses  contour  char's  of  constant  pH  at  5  selected 
depths  and  the  profile  is  generated  by  algorithm 
Contours  of  constant  pH  at  selected  depths  from  the  World  Ocean  Atlases 
of  Gorshkov  [  1 1 1  have  been  used  in  the  analysis  The  profile  increments  are 
0  1  pH  units  and  interpolation  is  required  /n-sftu  pH  values  are  used  in  the 
Global  Model  model  and  appropriate  corrections  have  been  made 
*  The  GEQSECS  reports  1 12)  contain  pH  profiles  obtained  from  circuits  of  the 

Atlantic,  Pacific  and  Indian  Oceans  These  data  have  been  used  together  with 
the  Russian  date  to  oevelop  the  constant-pH  contour  charts 
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Global  Model 


A=A,(MgS04)+A2(B(0H)3)+A3(rigC03) 

An=(S/35)0nP^fn/(Fn+»'^) 

a,  =  Fj -50x10^^^® 

r  ,/vT/70 

02=  0.1  K  F2==  0.9x10 

T/30 

83=  0.03  K  F3- 4.5x10 

Table  2:  Global  model  absorption  formula. 

VariabilUy  of  pH  is  clearly  the  major  limiting  factor  in  the  accuracy  of 
the  absorption  formula  In  Table  2,  the  pH  parameter  10^®^  has  been 
substituted  in  the  formula  of  Table  1. 

Five  K  points  have  been  found  to  give  a  reasonably  accurate  estimate  of  the 
K-profiles  The  points  are  numbet'ed  0-4  and  the  profile  is  fitted  by  the 
expression: 

K(0)=K{4)4  (Cq  ^  C,  D  ♦  Cj  ♦  C3  0^1  expl~(a  D)‘^) 

where  the  exponent  1  5  has  been  chosen  by  trial  to  obtain  'best*  results 
For  latitudes  less  than  60^,  the  depth  values  0=0,  0  5,  1,2,  and  4  km  are  used 
with  a=l/km  Since  water  depths  for  latitudes  greater  than  60°  are  generally 
shallower  and  most  of  the  variability  of  conceni  is  concentrated  nearer  the 
surface,  the  depth  values  0=0,  0  1 ,  0  3,  0  5,1  km  are  used  with  a=4/km. 

The  five  equations  for  n=0,1,2,3.4  to  be  solved  are  then  given  by 

*^0  *  *^1  *  *^4  &Xpi(aO^)‘ 

and  solution  for  the  coefficients  can  be  obtained  either  algebraically  or  by 
writing  tbs  equations  in  matrix  form  and  inverting 
The  K  contour  charts  for  the  lower  latitudes  are  shown  in  Appendix  A, 
Figures  1A-5A  The  high  latitudes  charts  ere  shown  in  Figures  6A-!  I A  and 
include  surface  values  for  both  summer  and  winter  conditions 
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K  Profiles 


Figure  6  shows  calculated  profiles  compared  to  typical  data  for  the  three 
regions.  Errors  are  seen  to  be  insignificant  in  comparison  to  the  expected 
overall  error  of  the  model.  The  North  Pacific  and  Eastern  Mediterranean  cases 
represent  the  expected  extremes  in  pH  effects. 

Both  K  and  temperature  p*"nfiles  are  used  for  integrating  absorption  along 
the  ray  paths.  The  temperature  profiles  are  generated  in  the  some  manner. 
Selinity  is  taken  as  constant  S=35  except  for  the  Mediterranean  case  where 
S=38  is  used.  Absorption  at  selected  frequencies  F(kHz)  at  ranges  R(km)  are 
calculated  by  integration  over  the  appropriate  ray  paths  for  the  convergence 
zctio,  bottom  reflection  and  surface  duct  modes.  Note  that  total  range  along 
tne  ray  paths  is  used  throughout. 

Expected  er»-ors  are  calculated  concurrently  using  AK=*0  05  for  the  CZ  and 
bottom-bouiico  modes  and  AK=±0.1  for  the  surface-duct  mode.  Values  for  the 
Thorp  equation  are  also  calculated  concurrently  so  as  to  minimize  relative 
errors.  The  values  A=A(Kmod)'A(Thorp)  are  the  dB  differences  between  the 
models  and  the  tA  values  indicate  the  degree  of  significance  of  the  K  model. 
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Convergence  Zone  Mode 


Atl  30"N 

Pac.  45'N 

E  Med. 

F  (kHz) 

R  (km)> 

74 

62 

39 

39 

K  mod. 

43.6 

29.8 

4.0 

Thorp 

40.4 

33.8 

oB 

21.1 

A 

3.2i1.1 

-4.0i09 

12.1±07 

5.5+0.6 

K  mod. 

mam 

25  0 

29  7 

23.0 

3.5 

Thorp 

34.2 

28.7 

17.9 

17.9 

A 

34tl.O 

-3.7i0  9 

1 1.8±0.6 

5.1+05 

K  mod 

32  0 

20.7 

26.2 

195 

3.0 

Thorp 

2b.7 

24  1 

15.0 

15.0 

A 

3.3±0.9 

-3.4i0.8 

1 1.2±0.6 

45+0.5 

K  mod 

26  7 

167 

16.3 

2.5 

Thorp 

23  8 

19  9 

12.5 

A 

2,9i0  8 

-3.3>0.7 

10  1+0  5 

38+0,4 

mod. 

21  5 

!5.0 

18  5 

13.2 

2  0 

rrv.,p 

192 

i6. 1 

10  1 

10,1 

A 

2.2i0.7 

-3  1±06 

0.4+04 

3.1+0  4 

Table  3:  CZ  two-way  path  absorption. 

Calculations  of  two-way  path  absorption  for  CZ  propagation  in  the  four 
regions  by  the  ray  integration  method  are  shown  in  Table  3. 

Values  for  the  N  Atlantic  indicate  that  the  K  model  is  roughly  3  dB  higher 
than  the  Thorp  formula  (Ai;3  dB)  in  the  3-4  kHz  range  This  is  due  due  to  fact 
that  K>1  below  the  axis  of  the  sound  channel  where  Thorp's  measurements 
were  made 

For  the  N  Pacific  case,  the  path  integrated  K<  1  and  As-4  dB  However,  the 
increase  in  K  with  depth  below  the  thermocline  causes  losses  to  be  greater 
than  for  the  sound-channel  mode  The  axial  value  Ks;0  5  would  give  a~-6  dB 
for  the  same  range  It  should  be  pointed  out  that,  in  some  areas  where  the  K 
profiles  do  not  fall  off  quite  so  rapidly  with  depth  within  the  theiTnocline, 
A-^0  dB  even  when  the  minimum  K  value  is  the  same  This  indicates  the  very 
important  effect  the  K  profile  and  the  hazards  of  trying  to  use  a  regional 
average  K  for  all  ray  paths 

For  the  Mediterranean  case,  K)l  over  the  entire  ray  path  and  A~i2  dB  The 
corresponding  value  for  the  Norwegian  Sea,  where  K  values  are  also  fairly 
uniform  but  not  as  great,  is  A;:5  dB  Note  that  CZ  conditions  are  probable  in 
both  regions  only  during  the  summer  months 
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Bottom  Bounce  Mode 


Rang.®  (km) 


Figure  7:  Bottom-bounce  ray  paths. 

Bottom-bounce  is  important  because  it  can  be  the  only  practical  mode  of 
operation  at  ranges  between  C2  zones  when  such  conditions  exist.  Figure  7 
shows  some  bottom-bounce  ray  paths  for  ranges  out  to  the  first  C2  in  the 
North  Atlantic  with  water  depth  4  km. 

Tables  4-7  show  calculations  of  the  Iv/o-way  path  absorption  for  the  four 
regions  considered.  Values  for  the  K  model,  Thorp  formula  end  the  difference 
A=Kmod-Thorp  are  shown  as  a  function  of  grazing  angle  and  ray-range  R(km). 
Note  that  K  model  values  are  approximately  equal  to  the  C7  values  corrected 
in  proportion  to  R,  This  is  exactly  true  for  the  Thorp  formula  since  the  loss 
IS  independent  of  ray  path. 

Bottom  grazing  angles  greater  than  about  5®  are  the  mein  concern  For  the 
N  Atlantic,  the  5®  K  model  values  ore  approximately  3.5  dS  higher  then  Thorp, 
Corresponding  A  values  ere  roughly  -45  dO  in  the  N  Pacific,  *9  dB  in  the  £. 
Mediterranean  and  *7  dB  in  tho  Norwegian  Sea  All  values  become  small  at 
high  grazing  angles  becousc  of  the  reduced  ronges 


2Alpha(clB)  BB  N.  At).  3.5kHz  4km  depth 


Theta" 

R(km) 

Kmod 

Thorp 

A(dB) 

5,3 

62 

32.1 

28.7 

7.2 

58 

30.0 

26.7 

3.3 

1 0.0 

53 

27.1 

24.1 

3.0 

140 

46 

23.5 

20.9 

2.7 

18,8 

39 

19  9 

n 

2.3 

24.0 

33 

16.9 

15.0 

1.9 

30.0 

27 

14.0 

m 

1.6 

36,0 

23 

11.7 

10.4 

1.3 

42.1 

mm 

9.5 

8.5 

1.1 

47.5 

8-1 

7.2 

0.9 

Tab's  4:  B3  two-way  path  absorption/N,  All. 


2Alph8(dB)  BB  N 

Par.  3.5kHz  4km  depth 

Theta" 

R(km) 

Kmod 

e 

o 

A(dB) 

5.4 

60 

27,4 

-4.5 

9.2 

53 

202 

omi 

-3,9 

143 

44 

17.0 

19,5 

37 

14.2 

1o9 

mt 

249 

32 

12  0 

143 

-2.3 

26 

99 

1  1  7 

-19 

36,5 

22 

62 

94 

-16 

419 

19 

mm 

0.5 

19 

47.1 

16 

6.0 

72 

■■■ 

517 

15 

5.5 

65 

■n 

Table  5  BB  two-way  oath  absorption/N  Pac. 
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2Alpha(dB)  BB  E.  Med.  3.5kHz  3km  depth 

Theta^ 

R(km) 

Kmod 

Thorp 

A(dB) 

5.3 

50 

31.9 

22.8 

9.1 

10.2 

41 

26.6 

18.9 

7.7 

16.0 

33 

BBi 

15.0 

6.3 

mm 

27 

17.6 

12.4 

5.3 

27.0 

23 

14.8 

10.4 

4.4 

32.9 

19 

8.5 

3.6 

38.4 

16 

10.2 

7.2 

3.1 

43.6 

15 

9.2 

6.5 

48.7 

12 

mm 

5.2 

53.2 

1 1 

6.5 

4.6 

1.9 

Table  6;  BB  two-way  path  absorption/E.  tlediterranean 


2Alpha(dB)  BB  70N10E  3.5kH2  3km  depth 

Theta* 

R(km) 

Kmod 

Thorp 

A(dB) 

5,3 

50 

29.7 

22.8 

6.0 

10.2 

41 

24.6 

10.9 

mm 

16.0 

33 

19  6 

15.0 

4.6 

mm 

27 

16.2 

12.4 

3.8 

27.0 

23 

13  6 

104 

32 

32.9 

19 

1 1.1 

8.5 

2.6 

38.4 

16 

94 

72 

2.2 

43.6 

15 

85 

65 

2.0 

487 

12 

6.8 

52 

1.6 

532 

11 

6.0 

4.6 

1.4 

Table  7  BB  two-way  path  absorpti on/Norwegian  Sea 
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•  CfQiinq  ongie,  Qeg 

Figure  8:  MGS  Bottom  loss  model 

Figure  8  shows  the  Marine  Geophysical  Survey  bottom  loss  curves  (or  the 
range  1-4  kHz  (131.  Province  numbers  characteristic  of  the  four  regions  are; 

N.  Atlantic,  *4.  N,  Pacific,  *6,  E,  Mediterranean,  *5,  Norwegian  Sea,  *3. 

The  MGS  data  analysis  was  based  on  the  Thorp  formula  and  the  question 
arises  about  the  possible  impact  of  the  K  model  A  good  indication  of  the 
approximate  magnitudes  can  be  seen  from  a  values  for  the  various  regions. 
Since  one-way  paths  were  involved  in  the  actual  experiments,  the  A  values  in 
Tables  4-7  must  be  divided  by  two 

The  results  indicate  that  K  model  effects  on  MGS  curves  are  significant  at 
the  lower  grazing  angles.  In  the  Norwegian  sea  and  E  Mediterranean  for 
example,  bottom  loss  should  be  reduced  by  some  3-4  dB  for  10®,  which  is 
roughly  equivalent  to  one  province  number  smaller  In  the  other  regions,  the 
corresponding  values  range  within  *2  dB  At  the  higher  grazing  angles,  the 
shorter  ranges  tend  to  make  the  effects  small  in  all  cases 
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Surface  Duct  Mode 


Frequency  (kHz) 


Figure  9:  Bismark  Sea  surface-duct  data  and  model. 


Figure  10  Gulf  of  Aden  surface-duct  data  and  model 
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Figures  9  and  10  show  attenuation  data  from  surface-duct  experiments  in 
the  Bismark  Sea  and  the  Gulf  of  Aden  [9]  compared  to  the  predictions  of  the  K 
model.  The  diffraction  loss  for  the  lowest  mode  is  calculated  for  a  bilinear 
gradient  with  g=.018/s  within  the  layer  and  a  large  negative  gradient  below. 
The  three  absorption  components  are  identified  and  the  top  curve  is  the  sum 
of  the  four  components. 

Since  the  model  accounts  for  the  experimental  data  above  the  diffraction 
cutoff  in  both  cases,  there  are  no  extra  loss  that  can  be  ascribed  to  surface 
reflection.  This  is  not  unexpected  because  the  sea  state  was  very  low  in  both 
experiments. 


2 Alpha  (d 

8)  SO 

Atl.  30®N 

Pac.45*N 

E.  Med. 

70N10E 

F  (kHz) 

R(kni)> 

25 

25 

25 

25 

4.0 

Kmod. 

20.1 

19.7 

mm 

mmm 

Thorp 

13.8 

13.8 

13.8 

13.8 

A 

6.3±0.9 

5.8i0.8 

8.0t0.9 

3.740.8 

3.5 

K  mod. 

18.1 

mam 

19.6 

15.1 

Thorp 

11.7 

1 1.7 

1 1.7 

117 

A 

6.4±0.8 

56±07 

7.81O.8 

3.4t0.7 

3.0 

Kmod. 

16.0 

WEom 

MBSm 

12.9 

Thorp 

9.8 

9.0 

9.8 

9.8 

A 

6.1i0.7 

5.2±07 

7.4tO.B 

25 

K  mod. 

■BH 

148 

10.7 

Thorp 

8.2 

8.2 

82 

8.2 

A 

5.4i0.6 

4.7i0.6 

6740.7 

2  640  6 

2.0 

K  mod. 

10.7 

10  6 

WEBm 

wmm 

Thorp 

6.6 

66 

,  66 

6.6 

A 

4  ltO.5 

40i0.5 

5.540.6 

2.1405 

Table  8;  Two-way  surface  duct  absorption. 

Table  8  shows  estimated  two-way  path  absorption  for  25  km  t'irget  range. 
Note  that  the  a  values  are  all  positive  in  this  case  because  pH  tends  to  be 
uniformly  high  near  the  surface. 

Since  surface  loss  must  also  be  considered,  the  question  then  arises  about 
the  reliability  of  the  prediction  models. 
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Frequency  (kHz) 


Figure  1 1;  Comparison  absorption  and  surface-loss  formula. 

The  empirical  surface-loss  formula  of  Schulkin  {14)  is  given  by: 

c4=l,6(Fh)’^^  d6/limiting-ray  cycle 

where  F  is  in  kHz  and  h  is  mean  waveheight  in  ft.  This  formula  was  derived 
from  analysis  of  AMOS  (15)  end  other  surface-duct  experiments  (161  by 
subtracting  the  S&M  absorption  formula  (li  from  the  measured  attenuation 
values,  dividing  the  result  by  the  skip  distance  of  the  limiting  ray  and  curve¬ 
fitting  the  results. 

Figure  1 1  compares  the  absorption  predictions  with  attenuation  estimates 
for  a  60m  surface  duct.  The  data  points  are  from  ref.  16.  To  estimate  the  K 
model  correction,  the  at,  and  S&M  formulae  are  added  as  shown  by  the  dashed 
lines  for  wave  heights  hrO  5,  1  and  2m  Subtracting  the  K  model  curve  from 
these  curves  then  gives  corrected  values  of  surface  loss.  Corrected  values 
are  therefore  smaller  by  an  amount  equal  to  the  difference  between  the  two 
absorption  formulae  in  this  example,  the  skip  distance  is  roughly  6.3  km  and 
the  difference  amounts  to  roughly  2  d8/bounce  at  3.5  kHz.  Since  the  K  model 
end  dashed  curves  are  nearly  parallel  in  this  range,  corrected  values  would 
become  negative  forh<  0  I5m,  approximately. 
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Figure  12:  Comparison  absorption  and  surface-loss  formula. 

Figure  12  compares  the  K-model  predictions  with  estimates  based  on  the 
empirical  surface-loss  formula  cujrently  used  in  the  GENERIC  Sonar  Model: 

as=-20  lQg,Q{0,3*0.7/lU(Fh/IO)^])  dB/limitlng-ray  cycle 

where  F  is  in  kHz  and  h  is  mean  waveheight  in  ft  This  formula  is  evidently  a 
modified  version  of  Beckmann-Spizzichino  theory  in  which  « (Fh)^  for 
small  parameter  values.  The  modification  causes  asymptotic  saturation  at 
ckjS  10  dB  for  high  parameter  values. 

The  dashed  lines  for  wave  heights  h=0  5,  1  and  2m  are  again  the  sum  of  the 
S&M  formula  and  ok,  for  the  limiting  ray  in  the  60m  duct  of  Fig  11  Where  the 
values  fall  below  the  K  model  curve,  correction  would  again  yield  negative 
values  of  surface  loss. 

It  IS  clear  that  both  surface-loss  models  tend  to  overestimate  surface  loss 
by  a  significant  amount.  However,  it  is  difficult  to  see  how  the  errors  can  be 
accurately  assessed  and  corrected  without  going  back  to  the  original  data  (if  . 
still  available !)  and  recalculating  surface  loss  using  the  present  absorption 
model 
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Conclusion 


This  study  has  been  limited  mainly  to  the  analysis  of  the  potential  impact 
of  the  K  model  on  active  sonar  systems  in  the  frequency  range  3-4  kHz  and  on 
the  surface  and  bottom-loss  models  used  in  performance  prediction.  The  four 
regions  examined  are  North  Atlantic,  North  Pacific,  Eastern  Mediterranean 
and  Norwegian  Sea. 

Comparison  of  K  model  predictions  with  those  based  on  the  Thorp  formula 
Indicates  that  two-way  path  loss  in  the  CZ  mode  can  vary  from  roughly  4  dB 
less  in  the  N.  Pacific  to  12  dB  more  in  the  E.  Mediterranean.  Exact  values  can 
depend  strongly  on  the  K  profile,  particularly  in  the  N.  Pacific  where  the  pH 
variations  are  greatest.  The  hazards  of  using  regional  average  values  of  K  for 
all  ray  paths  have  been  pointed  out. 

In  the  bottom-bounce  mode,  the  path  loss  is  roughly  equalivalent  to  the  CZ 
coefficient  times  operating  range.  Bottom  reflection  loss  must,  of  course,  be 
included.  The  MGS  bottom-loss  model  examined  here  was  based  on  analysis  of 
one-way  propagation  experiments  using  the  Thorp  formula.  In  the  Eastern 
Mediterranean,  the  corrected  bottom  loss  for  grazing  ang'es  near  10®  is 
roughly  equivalent  to  one  province  number  less,  the  amount  decreasing  with 
increasing  angle.  In  the  other  regions,  the  corresponding  values  are  generally 
smaller. 

The  K  model  predictions  for  two-way  path  loss  in  the  surface-duct  mode  at 
3-4  kHz  and  25  km  vary  from  3  to  8  dB  greater  than  Thorp  values  because  pH 
tends  to  be  highest  near  the  surface  The  higher  absorption  clearly  accounts 
for  a  large  part  of  the  observed  excess  attenuation  that  has  previously  been 
ascribed  to  surface  loss.  The  surface-loss  models  are  evidently  all  based  on 
analyses  of  surface-duct  experiments  using  pre-Thorp  absorption  models  and 
predictions  are  both  too  high  and  not  at  ell  consistent.  The  difficulties  in 
making  appropriale  corrections,  particularly  for  other  propagation  modes, 
are  obvious  and  a  major  reassessment  of  this  problem  is  therefore  indicated. 
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Figure  6A:  Surface  K  contours-summer 
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Figure  9A.  300m  K  contours 
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